‘e D 7) 
_ JUNE 74 
& B@ BOX 272 CALABASAS, CA 9313902 


IS 


VOL 2 NOG 


Ae 8 7°36. | 6 | 3 | 4 | | 1 
a ees ve | a | on 
ni We eS Se ee ee Cee eee ! 
ee2 | | tae 3. | : | 6 
ee 
| 4 | . | fi nell OD. 4 4 
| 2 | 2. g 
@ 3 & | eee 
om 9 3 
2 2 6 a 
“a 9 1 
i 
8B 


PC15-2 


Gane — re 


The POCKET CALCULATOR GAME is for two players, 
A and B. Players enter the board at the places marked, 
and proceed in turn, with the following move rule: 


Advance 1, 2, 3, 4, or 5 squares on the zigzag 
pattern. The score for each move is given by: 


NEEVINIVG: Ge 


where N is the number of squares advanced; C is the 
contents of the square landed on; and Co is the contents 
of the square the opponent last landed on. For the 
two initial moves, Co is zero. 


The game ends when one player overtakes or passes 
the other on the board, Ten points per square is added 
to the score of the player who passes the center square, 
outlined in red. 


EACH PLAYER KEEPS SCORE FOR HIS OPPONENT. Two 
markers, such as those included with the initial press 
run of this issue, facilitate keeping track of each J 
player's position on the board. The accompanying table 
of products of square roots is an aid to keeping score. 


The loser of one game starts first in the subsequent 
game, The sides of the board are alternated in 
successive games. 
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6 
2 1.4142 
3 1.7321 
4 2.0000 
5 2.2361 


‘N38 
tan! 8 


1.0000 | 1.4142 | 1.7321 | 2.0000 | 2.2361] 2.4495] 2.6458 | 2.8284 | 3.0000 


3.8730 | 4.4721 5.4772 | 5.9162 6.7082 


N-Series 4«8 
.60205999132796239042747778944898605353637976292422 
1. 386294 3611198906188 34464242916 3531361510002687 205 
1.5874010519681994747517056392723082603914933278998 
1.319507910772894259 37400197 12296401330334690131934 
1.219013654204475440911691002592560857 2774119358599 
1.148698354997035006798626946777927 58944 38508890978 
1..0139594797900291386901659996282 304258 363540227495 
1. 32581766 36680324650592392104284756311844406013064 


.90308998699194358564121668417 3479080304569644 38633 
2 .07944154.1679835928251696 364 3745297 042265004030808 
1.515716566510398082347 2598013064152 386812835429781 
1. 34590019263235613194283260375094.1904 3662470267775 
1.23114441334491628449939 306916774 31098761377611008 
1.02101212570719324976409517478306437 354108795327 37 
1.44644133224813518419996684.24758804.165254.145079177 


In the interest of completeness, the entries in the 
N-Series for 4 and 8 are given here. 
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Parkin 


Problem 15 called for the number of ways in which 
an 8 x 8 checkerboard could be cut into four congruent 
pieces, cutting along the lines of the checkerboard. The 
cover of issue No. 7 showed the 37 unique ways this could 
be done for a 6 x 6 board. 


Thomas R. Parkin, Vice President, Control Data 
Corporation, using an algorithm given below implemented 
on a CDC 6600, has extended the known results as follows: 


board size number of ways 
2 1 
4 5 
6 BiG 
8 782 
10 4ueko 


with an estimate for the 12 x 12 board of 7750000. 

Mr. Parkin estimates that the 12 x 12 case, using his 
algorithm on the CDC STAR, would take 5 to 10 minutes, 
and the 14 x 14 case from 30 to 60 hours. The following 
is reproduced from his letter. 
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It might be well, first. to show some comparison numbers to & 
support the seemingly large values for C{N}- We observe that 

each shape which cuts the square checkerboard into four 

congruent pieces is a polyomino of order N&. Therefore. we 

might suspect that the number CIN} would be less than the 

corresponding number of polyominoes of the corresponding 

number of squares. since. of course. there are poloyominoes 

of any given size {243 which cannot be used to divide the 

square array into four congruent pieces: i-e-.1 the 1 x Ib 


shape. Let's look at a comparative table: 
n ie 98 a 4 5 b 
N=2n a tn : 8 10 ihe. 
Ne ae ag by 100 444 
No/4SM Te ey 5 ite ag ae 
PLM} 1 § 2285 13079255 e 2.03x104° 2 §x1028 
CIN} ie = 5 ay 782 44240 ~ ?.75x105 


Thus we see that, in general, C{N} << PLN". 


On the other hand. one might argue that even though all the 
P£N°} cannot be used. one or more might be used several times; 


For example. the |_: - piece for N=4 can be used two different 


ways to yield two countably distinct divisions of the board into 
four congruent pieces. f{Incidentally. if I were setting the rules, 
I would exclude this case since the stated problem is "to cut 
the square checkerboard into four congruent pieces". and once 
these pieces have been cut; I maintain their origin is no longer 
germane’ however, this is not a problem except when N = 0 mod 4 
and it can be argued that the emphasis is on 'cut'. and not the 
shape of the cut piece, so I won't ieee the point-} Even so. 
an examination of the counts of P{N=} by symmetry type also 
shows that pieces possessing both 909 rotational symmetry at one 
point and 180° symmetry at another are vanishingly few, thus 

the conclusion, C{N} << P{N&} seems valid. 


There are a few observations about the problem which are 

necessary before one can understand the algorithm which I used 

for counting the possible cases. First of all. it is necessary 

to adopt a systematic procedure for enumerating the shapes. and, 

as a precursor to development of that systematic procedure. it & 
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is necessary to adopt some taxonometric scheme for classifying 
the shapes- Note that there are two classes of cuts which will 
divide the checkerboard into congruent pieces. and these may be 
classified as those which cut the square array into four pieces 
such that the edges of the pieces are 970° rotationally symmetric 
within the square: and the remaining cases where the square is 
divided into two rectangles by a line through the center, 
bisecting two sides. and then the rectangle is divided into two 
congruent pieces by a line which is 180° rotationally symmetric. 
Next; note that the pieces may be displayed either obverse or 
reverse, i-e-1 mirror images. or they may all be oriented in 
some canonical fashion. 


In my experience with writing programs to deal with two 
dimensional shapes, I have found that coding the shapes is 

quite tricky and difficult to work with’ however, dealing with 
the edges of shapes is sometimes simpler- Thus; in this problem 
I chose to consider the square array of 2N x 2@N squares in the 
plane to be a {2N+L} x {2N+]L} square array of lattice points 

in the plane and trace paths corresponding to the edges of the 
shapes of interest. 


The interior edges of the various pieces then provide for a 
taxonomy and allow a systematic procedure for generating shapes- 
Note that every piece in the square case involves an edge which 
traces a path from the outer perimeter of the square lattice 

to the center point. Furthermore. this path originates, says 
along the upper left border of the square down to the center 


line. and from no other place- {This is the canonical orientation. } 


It is these unique paths which connect the border with the center 
which we enumerate fin the square case} and similarly. with 
Suitable restrictions, for the rectangular case- 


We note. of course, that the origin of the paths in the square 
case ranges over only one-eighth of the border {because of 
rotation and reflection symmetry}. and over only one quarter 
of the border in the rectangular case- 


Now. of course, the algorithm is simple. Start at each 
appropriate border point and trace all possible unique paths 
to the appropriate center- I realize that this statement is 
descriptive but hardly constructive, so I will further detail 
the process. 


L.- Start at a corner; proceed down one side for across the topta 
and select the next point which has not been used as a 
starting point. Stop after using the center line point 
for less than or equal to one-quarter of the top}. After 
selecting the border pointi-mark all other border points 
as "disallowed". Enter the border point in a list at 
Level 1l.- 


2- Select the next lattice point inside the array. i-e-1 toward 
the center. and enter this point at Level 2- Note that 
there is only one possible choice for the second point on a 
path, given the first point on the border. 
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3. Mark the four for twot symmetrical points to the last & 
selected point as "disallowed". {The purpose of this 
disallowance is to preclude any point from appearing on 
more than one path, an obvious impossibility. thus. as each 
point on a path is selected. we simply check off the 
corresponding symmetrical points which are then no longer 
available for extensions of the current path-} 


4. Add to the list. at the current levels the coordinates of 
the 0, 1.1 21 or 3 possible points which could possibly be 
the extension of this path and mark them "unselected". 
Note that of the four points surrounding a given point. one 
of them is where the path came from. so that there are. at 
most. three possible extensions- Furthermore. it is possible 
to trace a path into a cul-de-sac such that no further 
extension of that path is possible and the path has not 
arrived at the center, hence the 0 possibility- 


5. If there are further eligible selections at this level, 
select the next unselected eligible extension point on this 
path at this leveli mark it "selected*., and enter it at the 
next level- If there are no further selections at this 
level, remove the "disallowed" exclusions for the selected 
point, and the selected point, at this level. Back up one 
level. If Level 1 has been reached. all done this border 
point’ if not. repeat this Step S. @ 


k. Match the selected point against the end point for pointst 
to see if this path is finished- If not. recursively 
repeat Steps 3, 41 5, and & until the path terminates or 
an exit at Step 5 occurs. 


?- Tally the terminated path fand note its coordinates. if 
you have time! } 


&- Backtrack one level and repeat at Step 5S. 


Using this algorithm. suitably modified in implementation to 

cover both the square and the rectangular cases. and again; 
suitably modified for the evenness or oddness of Ni yields the 
table of results, Table 1 attached. = 
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N 2N {Square} # {Rect} 
2 4 hea 2 ihe 
2.0 i O41 
3 b Tao 12 they 
240 q ed 
3.0 5 eth 
y 8 $0 242 140 
esl ib? 2.0 
3,40 146 3.0 
440 a? 4.0 
flea 
O.2 
5 “10 a0 Oso a 
2.0 B?74 Zab 
3,0 7538 3,0 
4.0 7224 440 
ceo usee 5.0 
ite 
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b Le 1.0 <1-8x105> 
: Ox 
bo0 <bxi0-> : 
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The Way to Learn Computing is to Compute 
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More Parkin 


Problem 30, The Web of Fibonacci (PC-10) called 
for the construction of a succession of triangles 
whose sides were the square roots of triplets from the 
Bioonacciwsequence (1, 1, 2; 3, 5, 8, 13, 21, 343...) 
The first point of the Web lies at (V2,0) and the 
second point at (V2,1); the problem was to find the 
coordinates of the 150th point. A solution by Thomas 
R. Parkin, Control Data Corporation, points out that one 
needs the value of Fj5,, where 


Fo Bij is) ody Hao 2: F3 = 3, and so on. 
Fisi is a 32-digit number: 
2609974810209 3884802012313146549. 


The 150th point lies on a circle of radius R = VFisi 2 


The pattern is formed of triangles whose central angle 
of interest is A = arctangent (Fi/Fra). 


"Thus it remains to compute the 148 A's, accumulate 
them, remove multiples of 27, compute X and Y from the 
remaining angle, and there we are. 


( 148 
"Thus A = >_arctangent (F,/Fisa) moa (27) 
ater 
A = 4.9844 radians 
R = en = ¥(26.09974810 x 102° 
R = 5.10879125646 x 1015 
X = R cosA = 1.37259283882 x 1029 


Y = R sin A = -4.92094879071 x 107° 


and we see that the point is just inside the fourth 
quadrant. All calculations were carried out in double 
precision Fortran on the CDC 6600 and are accurate to at 
least 12 significant digits for the final values of X 

and Y since the 96 bits allow for as much as 28 significant 
figures in all the intermediate calculations, including 

the trigonometric ones. 


"Finally we note that on approximately the scale 
of the cover of PC-10, the point in question would require 
a sheet of paper 200,000 million miles square to be plotted!" 
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Lake/fence 


PROBLEM 51 


A rectangular field is to be built with 500 feet of 


fencing next to a straight river. What are the dimensions 


of the field, if it is to have maximum area? 


The area, A, is X(500 - 2X). Setting the derivative, 

-2X + (500 - 2X), equal to zero, the length X is readily 
found to be 125 feet. This is not a computer problen, 
unless one seeks the required value by a bracketing process, 
to avoid using the calculus. 


If the straight river is replaced by a circular lake 
of radius one mile, the problem becomes more complex. An 
analytic solution is still possible, but difficult. see 
Figure D. ) 


H® = 52802 - (250 - x) 


2 


H® — 27815900 + 500x - x? 


t =tan7! ) (250 - x)/H 
Area A = Rectangle - sector OAMBO + triangle OAB 
Area A = X(500 - 2X) ~ 5280°(T) + H(250 - xX) 
and X is approximately 126.446... 
(A) Flowchart the logic of finding X to many places. 


(B) Write a program and run it, to find X to, say, 
15 significant digits. 
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500 - 2x 
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Comment 


Irwin Greenwald, Operating Systems Section, 
Xerox Corporation, writes: 


"In the 'Art of Computing 2! (PC11-4) you 
state 

'The test procedure should be logically 

independent of the program itself,'! 


"T beg to differ. While it is true that the 
test procedure should be concerned primarily with 
the function being computed, it will be less than 
satisfactory if it does not also account for the 
algorithm being used and idiosyncracies of the code 
itself. I am primarily concerned with the latter: 
all too often special cases, anomalous end conditions, 
etc., are introduced as a result of the manner of 
coding. Unfortunately, these are frequently known 
only to the coder, who may not be the person writing 
the test procedure. 


"As a case in point, consider the unnamed 
mathematician who expended considerable effort to 
test the precision of the JOSS LOG routine for values 
of the argument close to zero when, in fact, the 
agifficult precision problems occurred for values of 
the argument close to one. The 'tester' had no way 
of knowing that two different algorithms (and two 
different code sequences) were used depending on 
whether or not the argument was close to one. 
Obviously, a crucial part of testing this function 
was to select a range of test values around the 
‘erossover! points." 


One could hardly argue that it would not improve a 
program test if the person testing knew more about the 
program, But that wasn't the point. The point was that 
the test being devised should not repeat the logic of the 
program being tested. For a function subroutine like LOG, 
the test procedure should use some other algorithm, and 
for testing purposes that algorithm can be inefficient. 

Of course, the range of values should include end points 
and critical ranges. Perhaps a logically tight procedure 
would be to perform some extended calculation based on 
logarithms (e.g., the calculation of factorial 1000) that 
can be checked against the known result (PC2-10). 


ROBERT TEAGUE 


@ Speaking of Languages 


In issue No. 12, page 7, a set of three Fortran 
programs was presented as a test of how individual 
compilers interpret the ANSI standard. The programs 


10 


20 


are given again here: 


II. 


Til. 


DIMENSIGN L(3) 
EQUIVALENCE (J,L(2)) 


DIMENSIGN L(3) 


L(1) = 46 EQUIVALENCE (J,L(2)) 
DIMENSIGN L(3) L(2) 2 47 i eae 
EQUIVALENCE (J,L(2)) L(3) = 48 He = 9 
L(2) = 10 DG 10 J = 1,3 men) = 
Dg 10 J = 1,3 IeJ Dg 10 J = 1,3 
Ke Jd Tey er eee L(J) = J*4 + J/2 
CONTINUE 10 CONTINUE 10 CONTINUE 
WRITE (6,20) L(2) | write (6,20) L(2) WRITE (6,20) L 
FORMAT (17) | 20 FORMAT (17) 20 FORMAT (3(2X,17)) 
ST@P STOP STOP 
END END END 
Three questions were asked concerning each of 
them: 


(1) Will it compile? 
(2) Will it execute? 


(3) What output will be produced? 


In the outputs received, all three programs generally 
compiled and executed. This perhaps is a little surprising, 
since two of the programs, (2) and (3), reset the index 
value of the D@ statement. The EQUIVALENCE apparently 

very successfully hides this fact from the compilers. 


The outputs received are summarized in the table 


on the next page. 


Some interesting facts are evident fro 


m this exercise. 


The System 3 compiler apparently loads an index register 
and uses it as the D@ index without any reference to the 
variable J in storage. This is evidenced by the fact 

that although the value of J is changed to 10 during the 


second loop, and should therefore terminate 


the D@ since 


it then exceeds the test value of 3, the loop proceeds 
to a third iteration giving L(3) a value of 13. All the 
other compilers stopped after the second iteration, 


leaving L(3) with its original value of 3. 


The D&S 


‘Rt level compiler goes a step further, however, and tests 
for the D@ limit before the loop, since the initial value 


of L(2) is left unchanged, thus indicating a 
ation of the D@ is never made. 


second iter- 
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Compiler Output 
neal Gis ba Bact Pea 
1. System 3 3 3 4 3 13 
2. B6O700, Fortran 2.3 4 6 4 10 3 
3. IBM 1130, Eastern Michigan 4 6 4 10 3 
4, IBM 360/40, WATFIV 4 6 4 10 3 
5. IBM 360/40, 'F' Level D¢S 3 5 4 9 3 
6. CDC (without optimization) 4 6 Error Message 
7. CDC (with optimization) 1 7 Error Message 
8. CDC Run Fortran 3 5 5) ae 
9. NCR Full Fortran 4 6 5) 5) 5) 
10. Univac VM@S Background Fortran 4 6 5 S 3 
o 6 5) oes 


ee WATFOR Version 1, Level 3 


® 


The two runs on the CDC compiler certainly are of 
interest since the optimized version produces incorrect 
values while the unoptimized version produces values 
you would expect. However, it is notable that this is 
the only compiler (including WATFOR and WATFIV) that 
caught the fact that the index was being redefined. 


I would be interested in seeing additional outputs 
to these programs. Of particular interest would be 
outputs from 360/370 @S, and Honeywell compilers. 7 


Given below is the DATA DIVISION and PROCEDURE 
DIVISION for a CQBGL program. Add to the given code 
the IDENTIFICATIZN and ENVIRONMENT DIVISION entries 
necessary to run the program on your system, The 
question to be ascertained is this: what value will be 
printed for TOTAL? 


Readers have continued to work on the Change-Maker 
problem first presented in PC7-4 with an original solution 
presented in PC11-8. The first solution contained only 
9 Fortran statements. It now appears that the problem 
can be solved using only 8 Fortran statements and still 
conform to the ANSI standard requirements. Other short 
or novel approaches to this problem are welcomed. Send 
all material to "Speaking of Languages..." at Box 272, 
Calabasas, California 91302. 


DATA DIVISI€N. 
FILE SECTION. 
FD $UTPUT-FILE LABEL REC@RDS ARE @MITTED DATA 
RECORD IS SUTPUT-LINE. 
O01 MUTPUT-LINE. 
02 FILLER PICTURE X(5). 
02 TOTAL PICTURE Z(6).9(5). 
WORKING-STORAGE SECTIGN. 
77 SUM PICTURE 999. 
77 COUNT PICTURE 999. 
77 TEST PICTURE 999. 
77 IL PICTURE 999. 
PRECEDURE DIVISION. 
P-A, §$PEN ZUTPUT SUTPUT-FILE. 
MOVE SPACES TG QYUTPUT-LINE, 
MOVE ZERG? TS SUM, COUNT, TEST, 
P-B. PERF@RM P-C THRU P-D VARYING I FROM 1 BY 1 
UNTIL I = 8. GG TO P-E. 
P-C. ADD 1 T@ COUNT, ADD 1 TO SUM. 
IF I IS NOT EQUAL TS 4 OR 5 GS TH P-D. 
ADD I T@ SUM, ADD 1 T@ TEST, 
P-D. EXIT. 
P-E. COMPUTE TOTAL = SUM / (CQUNT ~ TEST). 
WRITE SUTPUT-LINE. CLOSE QUTPUT-FILE. 
STOP RUN. 


—_——Desk Caleulater Reviews ——— 


The Unisonic 739SQ (Unisonic, 16 West 25th Street, 
New York 10010) is a battery-operated, 8-digit pocket 
calculator with square root, reciprocals, % function, 
and one word of accumulating storage. On the rating 
scale published in issue No. 10, using the $68.88 price 
oe which the machine can be obtained, the rating is 
offdbo 


The Unitrex 800R (Unitrex of America, Inc., 11846 
East Washington Boulevard, Whittier, California 90606-- 
sold through Montgomery Ward and Company) is quite 
similar to the Uniscnic described above, but is a desk 
machine. It has square root, but no extra storage. 


Besides battery operation, an AC adapter can be obtained. 


At $59.88, its rating is 7.52. 


The ICP 537 ESR (ICP Ltd., New York 10001) is a 
scientific, battery-operated pocket calculator made in 
Taiwan. It is an 8-digit machine with floating point 
(but not scientific notation). There is one word of 
storage, and the following functions: log, in, 1/x, e%, 
square root, sine, cosine, tangent, arcsine, arccos, 
arctan, and the power function. If obtainable at $150, 
its rating is 9.33. 
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Log 15 1.1760912590556812420812890085 306222824 319 389827 285 
87 32351943817917812096350923661 355604110352943013 


Ln 15 2.7080502011022100659960045701487 1334417 30919120912 
6717 3647 34222511167 32809262667 3150374963290691170 


35 3.8729833462074168851792653997823996108 329217052915 
9082658757 376611348 30919 36979033519287 37685867 352 


{15 2 .4662120743304701014916113231545890427 354844866280 
5376017878741029337695292289746 3216298702664 34605 


{15 1.71877192758747877 7013521452044409157 1354589179517 
l 5367604292145116006836023906385989762028690950508 


415 1.472356700180346923712790000974024232796166407 5467 
7576 34682822148 393882655 3812353502808 303498650155 


15 1.3110194.23039749925204556407 052804 3307 320164 347835 
3539 31061269197 023347 2856344089 260863398747 237621 


yi5 1.0274505112667 269623447802620156116765316123602744 
58607 8000341600837407827 7593588060 376894 299644 342 


e 3269017 .3724721106393018550460917 213155057 385438200 
3420662956277 3242021 33274887913296987411229 


nt —- 28658145 .96938799845337882197 1660543332990696190337 
2849 3417989860167 1305896565088505771141500 


tan 15 1.504228163019072815032674997 3457803750007 114698620 
935407865158 309 330741106745767099 3225974688569116 


nN) 406561177535215237 397 279707 56704167 10103878906 32379 
7634290517698787 56383196170137717118109321745578199 
6250152587890625 


N-Series 
NO ee 


> 


